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Abstract 

An extended pion source, which can be temporarily created by a high en- 
ergy nuclear collision, will also absorb and distort the outgoing pions. We 
discuss how this effect alters the interferometric pattern of the two-pion mo- 
mentum correlation function. In particular, we show that the two-pion corre- 
lation function decreases rapidly when the opening angle between the pions 
increases. The opening- angle dependence should serve as a new means of 
obtaining information about the pion source in the analysis of experimental 
data. 

PACS numbers: 25.75.+r, 13.60.Le, 13.75.Gx 
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I. INTRODUCTION 



Hanbury-Brown-Twiss interferometry of identical particles emitted from a chaotic source 
has long been used to infer the spatial and temporal characteristics of that source 
||l|J^J^,|,|],P,^,H . Although originally conceived in astronomy to the end of extracting stellar 
radii the technique has found extensive application as a diagnostic tool in heavy- ion col- 
lisions ^,||J3J|J^. In most practical applications, the outgoing pions are assumed to suffer 
only weak final-state interactions. For instance, the interaction between the two detected 
pions is usually neglected, and a simple plane wave solution is used to calculate the two pion 
momentum correlation function, although in the case of charged pions, it is multiplied by 
the Gamow factor to correct for their long-range Coulomb force. 

In this paper, we discuss another type of final-state interaction, wherein each pion suffers 
in its exit from the source region. Namely, the interaction of the pions with the "source" 
itself. Since a pion source also absorbs and distorts pions, a source with finite extent in 
space and time will alter the outgoing pion waves in the source region. In particular, if the 
pion source is a droplet of "quark-gluon plasma," it may be regarded as a black-body of 
pions — and of all other hadrons — so that it casts its shadow on the paths of the outgoing 
pions emitted from the plasma surface. In a classical picture, this effect may be taken into 
account by introducing a strong correlation between the pion momentum and its emission 
points, so that the pions are emitted only in outward directions, an effect similar to that 
caused by collective flow [|I^,|rT[]. Pions may also interact in a dense pion medium via p or 
uj resonance formation, for example, and change the effective source size |]12[. One can also 



include baryon resonances explicitly in the pion source |13|. Yet there is a fundamentally 
more important effect which arises due to the quantum mechanical nature of interferometry: 
the distortion of the interference pattern generated by the pion-source-medium interaction. 

Here we study the pion correlation function in the presence of pion-source-medium in- 
teractions, as represented by a local optical potential, within the eikonal approximation. We 
calculate the distortion of the pion wave function explicitly, thereby allowing for a direct 
investigation of the quantum mechanical interference effects. We show that both the real 
and the imaginary parts of the optical potential give rise to non-trivial modifications: the 
former induces additional interference effects, whereas the latter leads to a suppression of 
the correlation as the opening angle between the pions increases. We demonstrate that 
the opening angle-averaged correlation function acquires a non-Gaussian shape due to these 
interactions, even if the source distribution is Gaussian. 



II. BASIC PRINCIPLES 

Before proceeding to include the final-state interactions, we first review the basic prin- 
ciples of two-particle interferometry. Such a description is easily found in the literature; we 
present it here to clarify its usual assumptions, as well as the manner in which we include 
pion-source-medium interactions. 

In essence, two particle interferometry results from a constructive — or destructive — 
interference of the two amplitudes for the emission of two identical bosons — or fermions — 
from two independent sources. Suppose that there are N such independent sources, localized 
at the space-time points Xi {i = 1, A^), and that the emitted particles do not interact with 
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each other after their creation. Let ip]i{xi) be the amphtude that a particle with momentum 
k is emitted from a source at x,. Then the pair amphtude \l/kik2(^i; ^i) that two particles 
are created by the source at Xj and Xj and are detected later with momenta ki and k2 is 
written in a product Ansatz as 



(1) 



We have taken out the overall phase factors associated with the arbitrary phase of each 
single particle amplitude, and since the detected particles are on shell, we can label the pair 
amplitude with just the three- momenta ki,k2. The first term is the amplitude that the 
particle of momentum ki is produced at Xi with the particle of momentum k2 produced at 
Xj. The second term is needed to ensure the proper symmetry of the total amplitude with 
respect to exchange of the particle coordinates — (+) for boson pairs and (— ) for fermion 
pairs. The probability of observing a particle with momentum k is given by 



A(k) 



1 

N 



(2) 



whereas the joint probability of finding two particles with momenta ki and k2 is 

2. 



P2(kl,k2) 



N{N- 1) 



^ ] ^kik2 (-^i) -^j , 



(3) 



Note that (...) indicates an average over the source distribution. The two-particle momen- 
tum correlation function is defined by 



C2(ki,k2) = 



P2(kl,k2) 
Pl(ki)Pi(k2) 



(4) 



With a few more assumptions, to be summarized below, the momentum correlation func- 
tion is related directly to the Fourier transform of the source distribution, Pq = J d'^xp{x)e^^^ . 
Namely, 



C2(ki,k2) = 1 ± \pq\^/\pq=0 



(5) 



where q = {Ei — i?2, ki — k2) is the relative four-momentum. Hence Eq. predicts that 

C2(ki,k2) ^ 2(0) as g^O (6) 



for boson (fermion) pairs, whereas C2 — >■ 1 as g — >■ cxd in both cases. 

There are three assumptions, besides the Ansatz in Eq. (P, involved in deriving the 
above from Eq. (HI). 



1. The relative phases of particles created at different locations are incoherent 
phase average is random. 



the 



2. There are no correlations in the source distribution. 
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3. The particles do not interact after their emission. 



The first assumption is the most important of all: it allows us to write the probability 
distribution in terms of an incoherent sum of the individual probabilities. That is, 

1 /■ 
^i(k) = j^^Y. WviA^i)?) = J d^xpiyx) |</?ki(a;)|^ , (7) 

P2(ki,k2) = ^^j^ _ ^^^ (El^kik2(a^»a^i)|^) 

= /<^'.,AJ(x,...)|*„„(x,...)|^ (8) 

where the averages over the point sources are replaced by integrals over smooth distribu- 
tions. All the interference terms containing the relative phase factor e**^''*"'^^-' drop out when 
averaged over the random phases 5i. Note that the two terms in the bracket of Eq. (|I|) have 
a definite relative phase, so that their interference survives the random phase averaging. 
The second assumption allows us to replace the two-particle source distribution ^2(3^1, 2^2) 
by the product of the single particle source distributions: 

p2{Xl,X2) = p{Xi)p{x2) . (9) 

This factorization does not hold in the presence of fine structure in the source, for example, 
in the multidroplet model of the mixed phase [0]. Lastly, by the third assumption, the 
single particle amplitude (p]s_{x) is given as a solution of the free Klein-Gordon or Dirac 
equation, as appropriate: 

M^) = Cue''"" ■ (10) 

Note that the normalization factors Ck cancel out in the correlation function of Eq. (^. 
Finally, these assumptions lead to the expression 

^ , X _ Jd^xi d^X2 p(xi)p(x2)|^kik2(a^i,a^2)P 

" /d%ip(xi)|y.k,(xi)|2/d%2p(x2)|^k2(a:2)r ^ ^ 



where 



^kik2 {xi, X2) = [ipk, (xi)(/?k2 {X2) ± V'k2 (a;i)^ki (^2)] . (12) 



III. INCLUSION OF PION SOURCE-MEDIUM INTERACTIONS 

We now turn our discussion to the modification of the correlation function due to the in- 
clusion of final-state interactions. In principle, there are two types of final-state interactions 
to be considered: the mutual interaction between the detected pions, and the interaction 
of the pions with the rest of the system. The former effect may be included by solving the 
Bethe-Salpeter equation or the Schrodinger equation for the two-body wave function 

Here we examine the latter effect, the medium modification [§|,|15|. We incorporate it by 
adding a self-energy term to the Klein-Gordon equation: 
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(□ + + 2mf/(x)) ipk{x) = 



(13) 



We have assumed that U{x) is a local optical potential. The heuristic derivation of Eq. (|TT]) 
in Sec. II presumes that the pions do not interact after their emission. However, Eq. (pUD 
still holds when a local one-body optical potential is included in the construction of (pk{x) 
as long as U{x) does not support bound states or create real particle pairs 

Many of the complications associated with the heavy-ion reaction are lumped in the 
potential U{x). For example, the pion source function p{x) appears explicitly in Eq. ([TI|), but 
the dynamical response of the source to the interactions of the outgoing pion with the source 
medium does not. This response is contained in U{x) — it is the effective single-channel 
optical potential. The final sum over the unobserved hadronic states implicit in Eq. (pUD 
can still be effected independently of the pion amplitudes in this picture. Furthermore, 
the initial state of the source function is also not identified, and for a given dynamical 
model description of the heavy-ion reaction, all allowed initial states must be averaged. We 
presume, for now, that this averaging can also be absorbed in U{x). 

The complexity of the optical potential is compounded by the time evolution of the 
source. For example, as the source expands, the relative momentum of the outgoing pion 
and of the source element with which it interacts changes in time and thereby alters the pion's 
interaction with the source. To illustrate this, consider an expanding source that consists 
merely of nucleons. The strength of the pion-nucleon interaction depends strongly on the 
relative momentum between the pion and the nucleon with which it interacts. Moreover, 
pion absorption by the source will depend on how two nucleons in the source are correlated 
in time. The complexity of U{x) can be reduced by constructing an explicit dynamical 
model of the heavy ion reaction, but simple physical insight into the effects of U{x) would 
be lost in such a task. Thus, we make the problem as simple as possible, in order to extract 
the essential physics that emerge as a consequence of the pion-source-medium interaction. 
To this end, we first choose to ignore the additional complications that arise from the time 
dependence of the source and consider merely a static pion interaction with the source 
medium. Note that Eq. ([I3|) may then be interpreted as describing the time-reverse of the 
usual scattering process: the pion returns from infinity with momentum — k and is inversely 
scattered by the optical potential U{x). For pions with energies of a few hundred MeV, 
ipk{x) is given approximately by 

M^) = Cke*'=--'Fr<^-'^(-'^'-') , (14) 

where {x, y, z) labels the pion production point. Note that z = k ■ r and z' parametrizes 
the straight line path along k from z to infinity. Equation ([l4|) is the eikonal approximation 
and agrees well with elastic pion-nucleus scattering data in the energy regime we consider 
||16| , when a realistic momentum-dependent potential is assumed. We will use the function 
(p-\i{x) in Eqs. (|lT]) and (0) to include the final-state pion-source-medium interaction in the 
pion momentum correlation function. 

Moreover, we choose a simple model for the pion-source-medium interaction — a complex 
constant multiplies the source density p: 

Tfi 

T^ir) = {an + iai)p{r) . (15) 
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In general, an and aj may depend on the momentum of the pion with respect to the medium 
and on the temperature of the medium as well. We shall assume that and cx/ are simply 
constant, but we will examine a wide range of parameter values in order to probe the 
sensitivity of the pions' momentum correlation to this new effect. We can then write the 
distortion factor fi'(k, r) as 

c/(k,r) = e-^("«+*"^)*(''''') , (16) 

where the thickness function t(k, r) = dz' p{h, z') is determined by the direction of the 
outgoing pion and its production point r. The impact parameter b is defined by r — z, where 
z = k ■ f k. The real part of the potential generates a phase shift, whereas the imaginary 
part represents a loss of flux. As we shall see below, both factors can significantly modify 
the correlation function. 

With the distorted single particle amplitude ipk{x) = e''^'^g{'k,r) as per Eq. (pIBD , the 
two-pion correlation function, Eq. (0), becomes 

C2(ki,k2) = l+|/(ki,k2)|2. (17) 

The second term, generated by the interference of the two distorted waves, is given by the 
square modulus of the overlap integral 

/(ki, k2) = ^ / d\p{r)g{k„ r)g{k2, r)e'^-^ , (18) 

where N = J d^rp{r)g{]<.,r). In the absence of pion-source-medium interactions, so that 
5f(k, r) = 1, the integral is reduced to the Fourier transform of the source distribution p(r). 
In this limit we thus recover the familiar result that the correlation function depends only 
on q = ki — k2. In general, however, /(ki, k2) depends on the relative angle of the two 
momenta ki and k2; hence, the correlation function also has this dependence. 

In Sec. Ill, we present numerical results for a finite, complex optical potential. However, 
a discussion of the black sphere limit (cx/ — oo, p(r) = poO{R — r)) serves as a useful 
illustration of the physics involved. Under these circumstances the outgoing pion amplitude 
attenuates to zero when its classical path transits the potential region. On the other hand, 
the source distribution p(r) is zero for r > R. We shall assume that the pions are emitted in 
a thin spherical shell of width 6 outside R and then compute the 6^0 limit. The distortion 
factor g{k, r) is determined purely by geometry: 




^(k, r) = I k ■ r + ^/ 1 - — ) ^(r - i?) . (19) 

Equation (p!8|) thus reduces to 

/(ki, k2) = ^ / dn 9ik ■ r) ^^(k2 ■ r) e"^-^ . (20) 

The resulting correlation function can be evaluated numerically for general 0. It can also 
be calculated analytically in the following two cases. For parallel pions, so that = 0,/ 
becomes 
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C2(k„k.2) = l + ^smM'^") . (22) 



SO that 



g2i?2 ^ 2 

In this case, the correlation function is simply a function of the dimensionless quantity qR, 
where q = |q| is the magnitude of the relative 3- momentum of the two pions. However, for 
pions emitted back-to-back, so that = vr, / vanishes. Thus, C2(ki, k2)|j^^^_j^^ = 1 for all q. 
The attenuation of the emission amplitude which penetrates the absorbing region destroys 
the interference entirely. 

The physical origin of the dependence can be readily understood in this limit and 
is illustrated in Figs. 1 and 2. For parallel pairs the deviation of the correlation function 
from unity is caused by the constructive interference of the two amplitudes corresponding 
to Figs, la and lb. However, for antiparallel pairs, the amplitude for the process in Fig. 2b, 
which would interfere with the amplitude in Fig. 2a, is completely suppressed by absorption 
in the source region; hence, in this case there is no enhancement. 

We can now proceed to evaluate Eq. (^Of ) for arbitrary and total pion momentum 
P = |ki + k2|. For given P and 0, there is a minimum q allowed kinematically; that is, 

/ 1 — cos (t)\ ^^"^ , ^ 

= ^ TT ! • 23 

\ 1 + cos ) 

As increases, gmin becomes non-negligible. At large opening angles, gmin is large, and only 
a small correlation between the pions remains — regardless of any pion-source-medium 
interactions. We therefore consider small opening angles, where gmin is small, so that the 
suppression of the correlation due to interactions with the source medium is easily detectable. 

In Fig. 3, we show the correlation function for several at two fixed values of PR in the 
above black sphere limit. As shown in Eq. (p^), the correlation function of parallel pairs 
goes to two at zero q. For comparison, the solid line shows the correlation function which 
results from ignoring final-state interactions. That is, (^(k, r) is set to unity in Eq. (|18|), and 
a uniform, spherical source of radius R is assumed, so that 



C^^'^iqR) = 1 + 



3 (sin qR — qR cos qR) 



(24) 



The two correlation functions are independent of and P, yet the scales with which they 
fall off in qR are different. In the black sphere case the correlation of parallel pairs at 
intermediate q is enhanced. This result can be understood physically as follows: when the 
source region is black, only pion pairs emitted from the front surface — as shown in Fig. 1 
— are detected; the pairs emitted from opposite sides as illustrated in Fig. 4 are suppressed. 
The effective source size is thus sma//er than it would be if all the pairs were sampled, leading 
to an enhanced correlation at intermediate q. When the opening angle between the pions 
becomes non-zero, the correlation function is suppressed relative to the parallel pion case, 
due to the final-state interactions. For PR = 10 (20), corresponding to P ~ 400 (800) MeV 
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for a source radius of i? ~ 5 fm, Fig. 3 shows the suppression of the correlation function as 
(f) increases from 0° (dashed hne) to 5° (crosses), 10° (squares), and 20° (pluses). 

Based on the above discussion, we hypothesize that these absorption effects should be 
quite general and thus present in fermion interferometry as well. 



IV. NUMERICAL RESULTS AND PHENOMENOLOGICAL CONSEQUENCES 

In this section we compute the correlation function ratio, Eq. (0), for the finite, complex 
optical potential of Eq. (plsl). For illustrative purposes, we choose simple forms for the source 
density p(r), a Gaussian p(r) = poexp(— r^/i?G^) and a uniform sphere p(r) = po9(Rs — r), 
with Rs = V2^Rq. The density parameters are chosen so that the two source functions have 
the same rms radii and central densities. Note that the integral / d?rp{r) differs in the two 
densities chosen. We equate the central densities of the two sources, rather than the number 
of particles / d?rp{r) they contain, as we believe the central density better characterizes the 
heavy ion collision. We presume that the fireball resulting from the heavy-ion collision is 
baryon-rich for simplicity, as is the case for heavy ion experiments in the nuclear stopping 



regime, such as those at the AGS [jT^,|T8[. We thus choose po = 0.17 fm~ . 

Given p(r), we need to estimate the constants aR and a/. In our treatment of U{x), aR 
and ai are proportional to the effective, forward vrA^ scattering amplitude in the fireball. 
The TiN amplitude in free space is well-known [|l^: due to 3-3 resonance formation, the 
imaginary part of the amplitude has a large peak near a pion laboratory energy of 200 MeV, 
around which the real part changes the sign. Below and above this energy region, the 
amplitude varies less and is of smaller magnitude. It involves many vrA^ resonances at higher 
energies. 

The construction of an effective ttN amplitude in nuclei from the free-space amplitude is 
an involved procedure, but the issues are well understood |20|. The main issues involved in 



its construction are as follows. 1) A resonance formation depends on the local relative mo- 
mentum of the 71 — N system, so that in the medium a minimally non-local optical potential 
is required. 2) Nuclear pion absorption involves more than one nucleon, so that the nucleon 
correlations in the medium become important. 3) Higher-order scattering contributions are 
now possible as well, and these also involve nucleon correlations. In view of the uncertainties 
in averaging over the initial source states and summing over the final states in constructing 
the potential in the fireball, we assume that these complications can be effectively included 
by varying the parameter values in our local potential. 

Since we use a simple model to illustrate the essential physics that emerge, we first 
choose a representative parameter set for an and u/ and vary it over a wide range. The 
representative set is based on 1) the free-space ttN scattering amplitude, taking into account 
2) nuclear pion absorption and 3) the non-observation of the initial and final states of the 
pion source, which acts to increase the pion's effective mean-free path pl|j2^ . We assume 



that the form of Eqs. (11) and (12) is preserved under the inclusion of this last effect, yet this 
need not be so. Finally, the representative set we choose is (ur, aj) = {±2, —2) fm^, which 
corresponds to a characteristic absorption length of (cxpo)"^ = 1.5 fm in nuclear matter. In 
order to explore a wider range of parameter values, we multiply this reference value by a scale 
factor c such that (aji, o"/) = (±2, — 2)c fm^ and vary c. To ensure that our examination is 
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not prejudiced by our choice of the imaginary to real ratio, we also examine the cases an = 
and CT/ = 0. 

In general, the correlation function depends on the kinematics of the pion pair — ki, k2 

— as well as on the parameters characterizing the source — R, a^, aj. For a spherically 
symmetric source, if final-state interactions are ignored, the correlation function depends 
only on q and R. However, if [/ 7^ 0, it depends on the opening angle and the total pair 
momentum P as well. In this section, we discuss these new dependences in turn. 

As discussed previously, we expect the correlation function to be suppressed at large 
opening angles as a consequence of including absorption in the source region. Our calcula- 
tions with a finite optical potential show that to be indeed the case. In Fig. 5 we plot the 
opening angle dependence of the correlation function at P = 400 MeV and g = 40 MeV for 
both a uniform source of Rs = 10 fm and a Gaussian source of Rg = 5 fm. The c = results 

— here there are no final-state interactions — do not depend on 0, yet there is always a 
significant suppression of C2 as increases for finite aj, in qualitative agreement with the 
black sphere results. We observe the same behavior if we fix aj and vary cr^, or vice versa, 
rather than varying them together. To wit, absorption by the source leads to suppression of 
the correlation at finite opening angles. 

If one averages the correlation function over the opening angle, the results turn out to 
be rather insensitive to P. We therefore fix P and investigate the sensitivity of the opening- 
angle-averaged C2 to the optical potential. The range of is restricted kinematically: 



< arccos 



P^-q^ 
P^ + q^ 



(25) 



Thus, for fixed P, as q increases, one includes an ever larger proportion of large opening 
angle pairs in C2, the correlation of which are suppressed. Angle-averaging thus results in 
a rather reduced correlation in the intermediate and large q range. At the lowest q, only 
parallel pairs are included, and the correlation approaches the Bose-Einstein limit, C2 = 2, 
as g — > 0, producing a non-Gaussian shape in the low q region. If one fits a standard 
Gaussian form to the correlation function, 

C2 = 1 + A exp [-g^Pe^l , (26) 



as is done with experimental data |T3], to extract the effective source size Re and the 
"coherence parameter" A, one finds a smaller effective source size Re < R and A < 1. The 
absorption of pions in the source region therefore leads to a small effective Re and A, as well 
as to a non-Gaussian shape of the correlation function in the low q region. 

In Fig. 6 we show the opening-angle-averaged correlation functions for both sources with 
several values of c. Note that P = 400 MeV and R = 10(5) fm for the uniform (Gaussian) 
source. For each source, we show the results for both signs of ct/j. Three curves are shown 
in each panel corresponding to c = (dashed lines), 1 (solid lines), and 10 (dot-dash lines). 
These results indicate a small but significant change in the shape of the correlation functions 
once the final-state interactions are turned on. Indeed, the medium- modified correlation 
function displays a non-Gaussian shape, especially at low q, even for a Gaussian source. 



Such non-Gaussian correlation functions have been seen |[T7|JTq.23]. For comparison, we 



show the E802 data for the vr correlations in 14. GeV/c ^^Si + ^^''Au collisions |T 
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together with our calculation in Fig. 7. The data can be fit to the form in Eq. (^) with 
Re = 3.42±0.26 fm and A = 0.65±0.07, yet the static Gaussian model we use also reproduces 
the correlation function if one assumes Rq ~ 5 fm and c 7^ 0. The solid line in the figure 
is calculated assuming P = 400 MeV and c = 10. The results for P = 600 MeV are almost 
identical. The agreement with data can be improved if we vary both the optical potential 
and Rg- However, our goal here is not a quantitative fit to the data — especially in light 
of our simple assumptions regarding the static source and its interaction — but, rather, a 
qualitative insight into how non-Gaussian shapes can be generated. Note that Fig. 7 also 
suggests that a study of the opening-angle dependence may be a better way to infer the 
extent of pion-source-medium interactions. 

If |(T/| is small, one may be able to detect the correlation effects induced by the real 
part of the optical potential. The extra phase introduced by the final-state interaction, 
exp{—iajit), modifies each of the pion wavefunctions and leads to an interference pattern in 
the correlation function. Such effects can be seen most clearly if a sharp edge is present in 
the source and the damping factor exp(— |cr/|t) is small. In Fig. 8, we show the 0-averaged 
correlation functions for a uniform spherical source with several values of and two values 
of 0"/. A large oscillation is observed in Fig. 8a when o"/ = and an is of roughly —1.5 fm^. 
The correlation function becomes as large as 4.5 at g = 100 MeV! Note, however, that we 
have not excluded bound-state forming parameters for < and small aj. The oscillations 
damp out quickly as ct/ grows large, as shown in Fig. 8b. 

At fixed P, q, and R, one can compare the values of C2 for different an and aj. Fig. 9a 
shows the variations of C2 with respect to ctr, for various fixed ai. In the upper panel, 
(T/ = 0, whereas aj = —2 and —4 fm^ — the solid and dashed lines, respectively — in the 
lower panel. The interference pattern is clearly seen at u/ = 0, yet only a small variation is 
observed for a more realistic value of aj. Note, however, the sudden jump in C2 when aji 
changes sign. Similarly, we study C2 as a function of aj for various in Fig. 9b. That is, 
we choose uj? = ±2 fm^ — the dashed and solid lines — and — the dot-dashed line. The 
interference induced by an enhances or suppresses the correlation for o"/? < or cr/j > if 
\ai\ is less than about 2 fm^. 

Even though our model is currently too simple to take its comparison with experiment 
seriously, the generation of modifications to the correlation function's shape by the pion- 
source-medium interaction is significant and robust. From our model the following simple 
picture emerges: the inclusion of pion-source-medium interactions leads to 1) a suppression 
of the correlation as the opening angle increases, 2) an angle-averaged correlation function 
which can deviate from a Gaussian shape at low q, as well as 3) possible additional interfer- 
ence effects from the real part of the potential. The last effect is suppressed by the inclusion 
of absorption, and when aj exceeds or is of the order of the value in normal nuclear matter 
the correlation becomes insensitive to the real part. If one fits the correlation data with a 
simple Gaussian form, Eq. (p6D, these effects will lead to a small apparent source radius and 
a small coherence parameter, A < 1. 

V. CONCLUSIONS 

Using a simple model for the pion-source-medium interaction, in the eikonal approxi- 
mation, we have demonstrated that an interesting momentum dependence in the two-pion 
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momentum correlation function is generated. Specifically, when absorptive interactions are 
finite, the pion correlation function decreases rapidly as the opening angle between the pi- 
ous increases. Our model is a simple one, yet this behavior should not be sensitive to the 
assumptions we have made. Pion-interferometry analyses have usually been performed as 
a function of only two of the two-pion momentum variables, implicitly assuming that the 
correlation function is independent of the opening angle. Our discovery of a strong opening- 
angle dependence thus suggests that analyses should be carried out explicitly in terms of 
the opening angle as well. Note, moreover, that the opening-angle-averaged pion correlation 
function exhibits a non-Gaussian shape when pion-source-medium interactions are included, 
even for a Gaussian source. This provides an alternative explanation for the coherence factor 
A, which is usually invoked to fit experimental data. 
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FIGURES 



FIG. 1. Emission of a parallel pair of pions from the source region. The amplitude in a) 
interferes constructively with the exchange amplitude in b). 



FIG. 2. Emission of an anti-parallel pair of pions from the source region. The amplitude in a) 
would interfere constructively with the exchange amplitude in b), but the latter is suppressed by 
absorption in the source region. 



FIG. 3. Two-particle correlation from a black sphere source of radius ii as a function of the 
relative momentum q = |ki — k2|. The total momentum P = |ki -|-k2| is 10/i? in a) and 20/i? in b). 
The solid lines indicate the results with no pion-source-medium interactions. The source region 
is assumed to be spherical and uniform. When the source medium is absorptive, the correlation 
function also depends on the opening angle (j) between the particles. Shown are the results for 
(p = Q° (dashed lines), 5" (crosses), 10° (squares), and 20° (pluses). For fixed P and (f), there is a 
minimum q, Eq. (^3|), allowed kinematically. 



FIG. 4. Emission of a parallel pair of particles from opposite sides of the source region. This 
process is suppressed compared to the emission in Fig. 1 if the source region is absorptive. 



FIG. 5. Two-pion correlation as a function of the opening angle between the particles. The 
relative momentum is fixed at g = 40 MeV, and the total momentum at P = 400 MeV. The optical 
potential is chosen with {(Jr,(Ti) = (±2, — 2)c fm^, and the results are shown here for a range of 
c. If the final-state interaction is ignored, so that c = 0, C2 does not depend on (j), as indicated 
by the solid lines. When absorption is present, the correlation decreases as (j) increases. Note that 
c = 1, 10 correspond to the dashed and the dot-dashed lines, respectively. Results for both an 
uniform sphere of radius Rs = 10 fm, in a) and c), and a Gaussian source of Rg = 5 fm, in b) and 
d), are shown with both signs of (Jr. 



FIG. 6. Two-pion correlation function averaged over the pions' opening angle, as a function of 
q. The total pion momentum is P = 400 MeV. Results for three different final-state interactions 
are shown, and, as in Fig. 5, c = 0, 1, and 10 correspond to the dashed, solid, and dot-dashed lines, 
respectively. 



FIG. 7. Comparison of the pion correlation function calculated using a static Gaussian source 
model [Rg = 5 fm) and the E802 data, corrected for acceptance and Coulomb effects, for the vr" 
correlation in 14. GeV/c ^^Si + ^^^Au collisions |18| with q. The solid line indicates the results 
using a = —20 — 20i fm^ (c = 10), whereas the dashed line represents the non-interacting (c = 0) 
results. In our calculation P = 400 MeV, but the results are insensitive to P. 
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FIG. 8. Opening-angle-averaged pion correlation function as a function of q for various cj/j and 
a uniform spherical source of Rg = 10 fm. P is fixed at 400 MeV, and in a) o"/ = 0, whereas in b) 
ai = —2 fm'^. The curves correspond to different values of aR] that is, the results with = —3, 
— 1.5, 0, and 3 fm^ are denoted by dashed, dot-dashed, solid, and dotted lines, respectively. A 
large oscillation is seen for gr ~ —1.5 fm^ if crj = 0, but very little sensitivity to ur is observed 
for more realistic aj. 

FIG. 9. Dependence of the opening- angle- averaged pion correlation function on a) gr and b) 
(T/. The calculation is performed with P = 400 MeV and g = 40 MeV for a uniform spherical 
source of Rg = 10 fm. In a), the upper panel shows the results for aj = 0, whereas the lower panel 
shows those for aj = —2 and —4 fm^ — the solid and dashed lines, respectively. In b), the results 
for aR = —2, 0, 2 fm^ are denoted by solid, dot-dashed, and dashed lines. The oscillations seen in 
Fig. 8 give rise to a large modification of the correlation at small \ai\. 



14 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 




(0 'Hd '}1^)^0 



(0 'Hd '}1^)^0 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



This figure "fig2-4.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 




+ exchange 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408005vl 



CD 
O 

:i 
o 

CO 

P! 

•i-H 

CO 
Ui 

:i 
O 



o 

V 

b 



/ 



/ / 

- / 
II / / 

I / 

I I ill I II I I I 



/ / 
/ 

/■ 2 
/■ II 

o 



o 

II 

o 

I I I I 



I I I I I I I 



o 

A 
b 



o 

II 
o 



I I I I I I I I I I I I I I 



\ / 

II ^ / 
^ I / 



/ 



/ 

/ • 

I II 
o 

/ 

/ 



I I I I I I I 



CL> 
CL> 

ft 



o 

• i-H 

d 



o 

V / 

K >i-l / . 

b II / / 

^ O 

* / 

/ 



V / 
/ 



/ 

/ / 



/ 



/ 



o 



I I I I II I I I I II I I I II I I I I I I I 



o 

A 

b 



o 
II 
o 



o / 



/ II / 
' I 



I I I I III I I I 



I I I I 



/ 

■/ 



o 

II 

- 

1 I I I I I n 



CO C\2 




^<(H'd'b)23> 



<(H'd'b)23> 



1 — I — I — r 




O lO O lO o 

q IN iq w o 

C\2 



<(H'd'b) 




^<(H'd'b)2o> 



^<(H'd'b)2o> 



